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We present an attempt at a critical overview of the current status of modeling for high-pr processes 
in nuclei. The paper covers several topics including coherence phenomena, in particular gluon 
shadowing and CGC; nuclear effects related to the restrictions imposed by energy conservation at 
large xl and xt ; space-time development of hadronization of highly virtual light and heavy partons 
and the related time scales; and the role of early production and subsequent attenuation of pre- 
hadrons in a dense medium. We identify several intriguing problems in the current paradigms for 
high-pT processes and propose solutions for some of them. 
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I. INTRODUCTION 



During the last decade of RHIC operation the physics of high-py reactions have attracted lot of attention, because 
these processes can serve as a probe for the early stage of the development of the quark-gluon matter created in heavy 
ion collisions. Intensive experimental and theoretical studies led to essential progress in our understanding in this 
field, especially for high-py processes on nuclei. Yet there is plenty of room for improvements and many experimental 
puzzles are waiting for solutions. This paper presents a brief critical overview of several key paradigms related to 
high-py physics. In particular, we challenge some of the current standard approaches, which have been settled by 
voting, rather than by real arguing. We also propose solutions for some of experimentally observed effects. 

In Sections 2-4 we start with the issue of the parton distribution functions (PDF) in nuclei, which are desperately 
needed for calculation of hard reactions at high energies. We present some of the theoretical tools for calculation of the 
gluon shadowing. Reviewing several popular gluon nuclear PDFs (nPDF), which resulted from global data analyses, 
we found some of them to be either ad hoc, or incorrect. The py dependent nPDFs are subject to the phenomenon 
of color glass condensate (CGC) controlled by the saturation scale. We demonstrate that the saturation scale in 
nuclei can be directly accessed by measuring py broadening in different reactions. Even more nuclear modifications 
affect the gluon distribution in the case of heavy ion collisions. We found a novel effect of mutual boosting for the 
saturation scales in the colliding nuclei, which can substantially enhance the effect of saturation in AA compared with 
pA collisions. 

In section 5 we study the restrictions imposed by energy conservation on hard reactions at large xy and xy, in the 
vicinity of the kinematic bounds. In view of the current controversy in interpretation of the suppression of high-py 
hadrons produced at forward rapidities in d-Au collisions at RHIC we propose a few tests, sensitive to the source of 
suppression, whether it is the deficit of energy, or coherence effects. Data available for these tests seem to favor the 
former mechanism. 

Section 6 presents a critical analysis of the current models for jet quenching observed at RHIC. We propose several 
tests for the popular scenario of energy loss, which is based on the unjustified assumption that hadronization is 
lasting a long time, and production of leading hadrons always occurs outside of the medium. We find that available 
data disfavor such a space-time structure of hadronization. Moreover, theoretically the production range of leading 
hadrons should be expected to be rather short because of the intensive vacuum gluon radiation ignited by the high-py 
scattering. 

In Section 7 we present a model based on our calculations for the pre-hadron production length. We demonstrate 
that the produced pre-hadron dipole quickly expands its size, especially if it contains an open heavy flavor. This leads 
to a strong absorption in a dense medium created in heavy ion collisions. In the limit of short mean free-path of such 
dipoles the nuclear suppression factor Raa for high-py hadrons can be predicted in a parameter free way. RHIC data 
for light and heavy flavored hadrons at high py seem to support existence of such a regime. 

The main observations are summarized in section 8. 



Bound nucleons are well separated in the nuclear rest frame. They still do not overlap in the nucleus Lorentz boosted 
to the infinite-momentum frame. Indeed, although the Lorentz boosted nucleus looks like a pancake, not only the 
nucleon spacing, but also the nucleons themselves are subject to Lorentz contraction. However, partons carrying a 
small fraction x of the nucleon momentum are contracted x times less. As a result the parton clouds originated from 
different nucleons overlap in the longitudinal direction at small x <C 1 [l] , as is illustrated in Fig. [TJ 

This leads to a dense packing of the radiated gluons in the phase space. However, according to the Landau- 
Pomeranchuk (LP) principle Q multiple interactions do not generate multiple radiation of gluon spectra, if the 
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FIG. 1: Left: Bjorken x (vertical axis) dependence of gluon clouds spread over longitudinal coordinate (horizontal axis) in a 
high-energy nucleus. The Lorentz factor decreases downwards small x, and the clouds start overlapping. Right: Fusion and 
interaction with the target of gluons originating from different bound nucleons. The amplitude squared gives the triple-Pomeron 
graph. 



coherence time of radiation is large, l c ^> Ra- The radiation process does not resolve between single and multiple 
interactions, only the total momentum transfer matters. In other words, when the interaction becomes sufficiently 
strong, it starts screening itself. In particular, radiation of soft gluons, which has large cross section, can be strongly 
shadowed. In terms of the Fock state decomposition, a gluon in a given Fock state can be radiated only once. 

Therefore, the spectrum radiated with small transverse momenta, k\ <C (k^) from multiple interactions must 
saturate when the phase space of radiated gluons is densely packed. Only at sufficiently large transverse momentum 
of gluons, kx > (kr), where the phase space becomes dilute, multiple interactions start contributing to the multiplicity 
of gluons, increasing the range of kr- Eventually, one arrives at the Bcthc-Hcitler regime of radiation when each of the 
multiple interactions equally contributes to the radiation spectrum. The transverse gluon momentum characterizing 
the transition scale between the two regimes is called saturation momentum and is defined below. This effect is known 
as the color glass condensate (CGC) @-|5|. 



A. CGC: how to measure the saturation scale 



The rise of total cross sections with energy, discovered back in 1973, was the first manifestation of the increasing 
population of partons at small x. As the parton densityincreases, the inverse process of parton fusion becomes 
important, and eventually the parton density is expected |6| to saturate. 

In the nuclear rest frame, the same phenomenon looks like Glauber shadowing and color filtering for a dipole 
(quark-antiquark, or glue-glue) propagating through nuclear matter [7J • The partial elastic dipole- nucleus amplitude 
at impact parameter b reads [8| 

/£(&) = l-e-H>-*)^), (1) 

where tt and E are the transverse separation and energy of the dipole respectively; 1^(6) = J^^dz pa(o, z) is the 
nuclear thickness function, integral of nuclear density along the trajectory of the projectile at impact parameter b. We 
assume here that the Lorentz dilated length (coherence length) of the dipole size fluctuations is much longer than the 
nucleus. The qq dipole cross section on a nucleon should vanish at small transverse separation o~q q {rT — > 0,E) oc 
Q. The energy dependence is discussed below. For large Ta(6) only the small- tt part of the cross section contributes, 
so one can use the r\ approximation, 

a» q (r T ,E)\ rT ^ vC q (E,r T )r 2 T , (2) 

where C q (E ) rr) is logarithmically divergent at small tt [1]. I n this limit the factor C q can be related to the gluon 
distribution [9(, 

C q (E, r T ) = a s (l/r 2 T ) xG(x, 1/4), (3) 



where 1/x = Im^E ' r\. 
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The quark saturation momentum Q q A is usually introduced as, 



/£(&) = 1-exp 



r$Q 2 qA (b,E) 



so comparing with Eqs. ([T]) and © one gets, 



f qA (b,E) 



2C g (E,r T = l/Q q A)T A (b). 



(4) 



(5) 



Here we fixed the dipole separation at the typical value ~ l/Q qA relying on the weak, logarithmic, tt dependence 
of C(E,r T ). 

The energy dependence of the factor C(E) fitted to DIS and photoproduction data is depicted in the left panel of 
Fig. [2] It steeply rises with energy up to very high values. We remind that because Q qA increase with energy, the 
mean value of {tt), which we fixed above, should be reconsidered. 
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FIG. 2: Left: Factor C(E) defined in ((2}, which controls the saturation scale in nuclei through Eqs. l[6])-([7|). The dashed curve 
shows the dependence on quark energy in the leading order approximation. The solid curves show the modified broadening 
factor C q (E) — R g (E) C q (E) darned by gluon shadowing, which depends on nuclear thickness Ta propagated by the quark. The 
curves from bottom to top correspond to Ta = 1.5, 1.0, 0.5 fhi -2 . Right: Broadening in Drell-Yan reaction on different nuclei 
as measured in the E772 (closed squares) [Jj| and E866 (open squares) [Jj|] experiments, respectively. Broadening for J/\t and 
T [l2l [Til ] is shown by circles and triangles, respectively. The dashed and solid curves correspond to the predictions without 
and with the corrections for gluon shadowing, respectively. 



The same factor C{E), as in Eq. ([5]), controls broadening of the transverse momentum of a single parton propagating 
through a nucleus [lja, [Tj| , 

= 2T A {b)C q (E,r T = 0). (6) 

Thus, we arrive at a divergent result, since C{E,tt) ^ at r-r — > Q. This is not a surprise, as the mean 

transverse momentum s quar ed (jpj) is ultraviolet divergent. Moreover, this divergency is not canceled in broadening 
Ap|, = (pj) A — {Pt)n I id fTT| . To settle the problem, one should fix tt at a characteristic value like in ([5]) at 
r% ~ 1/Ap|. 

Thus, broadening and the saturation momentum are equal: 

Q 2 qA (b,E)=ApUb,E), (7) 

so one has a direct access to the saturation scale by measuring broadening. 

The experiments at HERA provided detailed information about the dipole cross section a qq (rT,x) as function of 
x and the dipole size. Therefore the factor C(E,rr) in ^ and broadening can be well predicted. The results are 



A P * T = T A (b) V r 2 T <(r T ) 
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depicted by the bottom dashed curve in Fig. [5] (right) in comparison with broadening in Drell-Yan reaction measured 
in the E772 and E866 experiments [HI, IHJj- Notice that this is broadening or the saturation scale for quarks. For 
gluons it is expected (perturbatively) to be 9/4 times larger. That should be compared with broadening of heavy 
quarkonia. Comparison of the predicted broadening shown in Fig.[2]by the two dashed curves, with data for J/ty and 
T [12|, [lj] also demonstrates good agreement. 

Notice that to calculate nuclear broadening for heavy quarkonium production one does not need to know its 
mechanism provided that the coherence length of quarkonium production is short, i.e., 

V!l «Ra- (8) 



m N M% 



The data for broadening of T depicted in Fig. [5] satisfy well this condition, only the data for J/ty production arc 
somewhat out of this kinematic domain. 



III. GLUON SHADOWING 



As was explained above and is illustrated in Fig. [TJ the overlap of gluons originated from different nucleons leads to 
their fusion and to a reduction of the gluon density at small x. This is how gluon shadowing looks like in the infinite 
momentum frame of the nucleus. In the nuclear rest frame the same phenomenon is interpreted quite differently It is 
related to higher Fock components in the incoming hadron, containing gluons. Multiple interactions of those gluons 
give rise to gluon shadowing. This interpretation has the advantage of being more intuitive and directly linked to the 
optical analogy of shadowing. 



A. Theoretical expectations 



Unfortunately, no satisfactory theoretical description of gluon shadowing, which would work in all kinematic regimes, 
has been developed so far. Each of the existing approaches has limitations, which makes a direct comparison with 
data difficult. 

The most rigorous quantum-mechanical treatment of gluon shadowing is based on the path-integral technique 
fl7l HH . It includes both the phase shifts between different production points and a proper treatment of the dipole 
attenuation. Therefore this approach should correctly reproduce the onset of shadowing at medium-small x where 
the magnitude of shadowing is small. Actually, shadowing is expected to remain weak even at smaller x because of 
the specific nonperturbative effects suppressing gluon radiation [17| . Indeed, high-statistics data for single diffraction, 
pp — > pX, show that gluon radiation, which controls the large invariant mass tail, da s dj dM\ oc l/M^, is an order 
of magnitude smaller that one could expect inpQCD. This problem has been known in the Regge phenomenology 
as smallness of the triple- Pomeron coupling [19]. The only way to suppress gluon radiation is to reduce the mean 
quark-glue, or glue-glue separation. In order to explain data for diffraction one has to assume this separation to be as 
small as ro ~ 0.3 fm [llj ■ Indeed, lattice calculations lead to a glue-glue correlation radius of this order of magnitude 



|20j . There are many other experimental evidences supporting the existence of such a semihard scale in hadrons [21 

Thus, the magnitude of gluon shadowing evaluated in [l7j is rather small, as is confirmed by calculations performed 
in [22| and depicted in the left panel of Fig. [3] 

As we mentioned above, this is the most accurate method, which is valid in all regimes of gluon radiation, from 
incoherent to fully coherent. Nevertheless, this is still the lowest order calculation, which might be a reasonable 
approximation only for light nuclei, or for the onset of shadowing. Contribution of higher Fock components is still a 
challenge. This problem has been solved so far only in the unrealistic limit of long coherence lengths for all radiated 
gluons, in the form known as the Balitsky-Kov cheg ov equation (BK) [23|, [24j. Numerical solution of this equation is 
quite complicated and includes lot of modeling [25] ]. A much simpler equation, which only employs a modeled shape 
of the saturated gluon distribution, was derived in (26[. It leads to a gluon distribution in nuclei, which satisfies the 
unitarity bound I27T| , and the results are quite similar to the numerical solutions [25[ of the BK equation. The new 
equation reads [26|, 

Rl 



= ' " (l + R g n ni+n eff ) (9) 
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FIG. 3: Left: The nuclear ratio R g for gluons as function of path length in nuclear matter, calculated in [221 ] at Q = 4GeV 
for several fixed values of x. Right: R g as function of nuclear thickness Ta ~ 0.16 fm -3 x L, calculated at the semihard scale 
Qo — 1.7GeV 2 and x = 10 -4 . The dashed curve presents gluon shadowing corresponding to hadron-nucleus collisions. The 
solid curve includes the boosting effects specific for central nucleus-nucleus collisions (Ta = Tb). 



where 



n (E,b) 



9C(E) 
2<4v(£) 



n eff (E,b) = -C(E)r 2 T A (b). 



(10) 



The energy dependent factor C(E) was introduced above in ©-{HI- The mean size of a gluonic dipole, ro sa 0.3 fm, as 
was already mentioned, is dictated by data [I3,[Hi- We rely on the saturated shape of the dipole-nucleon cross section 
with the saturation scale Q qN (E) = 0.19 GeV x (E/lGeV) 014 , fitted to DIS data [HI HI]. The solution of equation 
© at x = 10~ 4 and Q% = 1.7 GeV 2 [2f| is depicted by dashed curve in the right panel of Fig. [3] Compared to the 
results depicted in the left panel we see that Eq. © predicts a somewhat stronger shadowing than the path-integral 
method (Q 2 dependence is rather weak). This is not a surprise, the left figure was calculated accounting for phase 
shifts, which are not yet small even at x = 10~ 4 . Indeed the x-dependence of shadowing shown in the left panel does 
not seem to saturate at this value of x. This means that in the saturated regime with very long coherence lengths for 
all gluons, which was assumed for the equation ([9"]). shadowing should be stronger. 

Notice that Eq. ([§]) does not contain any hard scale, but only a semi- hard one controlled by ro [2l|. Therefore 
its solution should be treated as the starting gluon distribution at the semi- hard scale Q\ ~ 4/ r o- Then shadowing 
should be DGLAP evolved up to an appropriate hard scale. 



B. Gluon shadowing from data 



Since gluons dominate at small x, the nuclear modification of the gluon distribution is desperately needed, if one 
wants to understand nuclear effects in hard reactions at high energies. There is wide spread believe in the community 
that gluon distributions in nuclei are known from available global analyses of data. Unfortunately, we should admit 
that most of the nuclear gluon PDFs, which are currently used in calculations and data analyses, are either an 
educated guess having little too do with data, or just plain wrong. 

As an example, data for J/ty production from the PHENIX experiment were analyzed [29[ using two popular gluon 
PDFs, which according to the authors of [29| come from data: "The modified nuclear PDFs from EKS and NDSG 
are constrained from other experimental measurements such as deep inelastic scattering from various nuclear targets 
and the resulting Fz(A) structure functions." Well, let us check whether such a claim is really justified. 

First, the nDSg plotted in the left panel of Fig. SI which presents the result of the global analysis of DIS data 
performed in [30] . The gluon ratio R g which comes from the data analysis is depicted by the upper dashed curve, 
while the bottom dashed curve, labelled by nDSg, is an ad hoc curve, which actually contradicts the DIS data. This is 




FIG. 4: Left: Nuclear ratios at Q 2 = 10 GeV 2 as function of x fitted in 30] to DIS data. The upper dashed curve presents 
the results of the fit for the gluon ratio R g . The bottom dashed curve (nDSg) is an ad hoc trial function (see text). Ri ght : 
The gluon ratio R g provided by the EPS08 code at Q 2 — 1.69 GeV 2 . The shaded area is forbidden by the unitarity bound [23] 
shown by dashed curve. 



what the authors of [3(| say about it: "We provide the result in set called nDSg, constrained to satisfy Rg u — 0.75 at 
x N = 0.001 and Q 2 = 5. The x 2 value of this analysis is around 550, considerably larger than the unconstrained fit, 
and should be considered only as a mean to study variations on, mainly, the gluon nuclear distribution." Strangely, 
the authors of [29| ignored the real solution for R g , but picked up the fake curve, which provides a several times 
overestimated gluon shadowing. 

The magnitude of gluon shadowing in EKS98 is similar to nDSg, and does not come from data either. The gluon 
shadowing at small x also was fixed at an ad hoc large value, R a = Rf 2 , what resulted in an overestimated gluon 
shadowing and antishadowing. This is what the authors of [3l| say: "We point out, however, that even though we 
obtain a very good agreement with the NMC data [3] and with the analysis of Ref. [5], we can confirm our initial 
assumption of gluon shadowing at small values of x only on fairly qualitative grounds (stability of the evolution), 
rather than through a direct comparison with the data." 

The recent new global analyse s 13311 tried to enhance the input from data by including BRAHMS results for hadron 
production at forward rapidities [34| . This was a risky addition, because the mechanisms of suppression are still under 
debate (see Sect. [Vj, and no consensus has been reached so far. This was probably the reason why this fit led to 
such astonishing results. The gluon shadowing was found incredibly strong, stronger than is allowed by the unitarity 
bound (27| . as is demonstrated in the right panel of Fig. |4j Apparently, the problem is caused by misinterpretation 
of the suppression observed in [34| relating it to gluon shadowing. 

Indeed, the later analysis EPS09 35], which wisely excluded the BRAHMS data, led to the PDFs, which failed to 
reproduce the suppression observed in the BRAHMS experiment at 77 = 3.2. It disagrees even more with STAR data 
at r\ = 4 [36|] . However, in the new analysis EPS09 the BRAHMS data was replaced by a new source of information 
about gluon nPDF, the Cronin effect in neutral pion production at the mid rapidity at RHIC, which created a new 
problem. The mechanism of the Cronin enhancement was quite misinterpreted in 35] . This effect is a manifestation of 
the CGC phenomenon, which is a modification of the transverse momentum distribution of partons in nuclei increasing 
their mean transverse momentum up to the saturation scale. This also can be understood as broadening of partons 
propagating through a nucleus, as we discussed in Sect. [TTJ 

Such a modification of the p^-distribution in nuclei has been always considered as a mechanism for the Cronin 
effect. This mechanism explains well data from fixed target experiments [37l l38l and correctly predicted [37j the 
Cronin effect for pions measured at RHIC [39| . These data were interpreted in [35] in the collinear approximation, 
but the Cronin exhancement was attributed to the nuclear modifications in nPDFs. This mechanism can be disproved 
comparing it with data from fixed target experiments, where the Cronin enhancement is huge and cannot be explained 
by nPDFs. So we conclude that gluon nPDF in EPS09 are not trustable, since the main source of information about 
nuclear modification of the gluon distribution is based on the incorrect dynamics for the Cronin effect. 
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We should also comment on the global analysis HKN07 which was performed at the NLO level. Among different 
kinds of modeling in the current theoretical estimates of gluon shadowing, some are more debatable, some less. In 
particular, our knowledge of the coherence time of gluon radiation, which controls the onset of gluon shadowing, is 
pretty solid. Shadowing is possible only if the coherence time becomes longer than the mean internucleon spacing, 
about 2fm. It is clear that a \q~qg) fluctuation of a virtual photon is heavier than a \qq) one, therefore the coherence 
time for gluon shadowing must be shorter than for quark shadowing. This is a very solid, non-debatable statement 
dictated by kinematics. Thus the onset of gluo n shadowing must occur at smaller x than for quarks. The latter is 
well known both from theory [4l[ and data |42j to occur at x ~ 0.07. 

Such a restriction should be used as a physical input when the initial x-dependences of the PDFs at the starting 
scale are shaped. This apparently has not been done in the analysis [40] , which led to an unphysical x-dependence of 
gluon shadowing, which starts at x = 0.2 and is rather strong already at x = 0.1, where it is excluded by the above 
arguments. 

IV. BOOSTING THE SATURATION SCALE IN AA COLLISIONS 

Due to broadening a nuclear target probes the parton distribution in the beam hadron with a higher resolution. 
Therefore, the effective scale Q 2 for the beam PDF drifts to a higher value Q 2 + Q 2 ^. At first glance this seems 
to contradict casuality, indeed, how can the primordial parton distribution in the hadron depend on the interaction 
which happens later? However, there is nothing wrong. The interaction performs a special selection of Fock states in 
the incoming hadron. The same phenomenon happens when one is measuring the proton parton distribution in DIS. 
The proton PDF " knows" in advance about the virtuality of the photon which it is going to interact with. 

The shift in the scale also can be interpreted as a manifestation of the Landau-Pomeranchuk principle Q : at long 
coherence times gluon radiation (which causes the DGLAP evolution) does not depend on the details of multiple 
interactions, but correlates only with the total momentum transfer, q + Apx, which after squaring and averaging over 
angles results in Q 2 + Ap\. 

As far as the PDF of the projectile proton has a harder scale in pA collisions than in pp, the ratio of parton 
distributions should fall below one at forward and rise above one at backward rapidities. This may look like a 
breakdown of fc^-factorization, however, it is a higher twist effect. 

Examples of pA to pp ratios Ra(x, Q 2 ) calculated with MSTW2008 [43| are shown in Fig. [5] for gluon distributions 
in a hard reaction (high-p^, heavy flavor production, etc.). We see that the shift in the hard scale caused by saturation 




FIG. 5: Left panel: Nucleus to proton ratio of the gluon distribution functions in a reaction characterized by the saturation 
momentum Q\a = 2 GeV 2 and hard scale Q 2 = 2, 3, 5, 10 GeV 2 . Right panel (top): pA collision in which the colliding proton 
is excited by multiple interactions up to the saturated scale Q 2 s a, which leads to an increased multiplicity of soft gluons in 
the incoming proton. Right panel (bottom): nuclear collision in which participating nucleons on both sides are boosted to the 
saturation scales, Q\a i n t ne nucleus B, and Q 2 sB in the nucleus A. As a result, the low- a; gluon population is enriched in both 
nuclei. 

in the nucleus leads to a sizable suppression in the projectile parton distribution at large x — > 1 and enhancement at 
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small x <C 1. We also observe that the magnitude of nuclear modification quickly decreases with Q 2 confirming that 
this is a high twist effect. 

Notice that in pA collisions the modification of the PDFs of the beam and target are not symmetric. Namely, the 
scale of the PDF of the beam proton gets a shift, Q 2 => Q 2 e ff = Q 2 + Q 2 s a-> while the PDFs of the bound nucleons, 
which do not undergo multiple interactions, remain the same as in pp collisions. This is illustrated pictorially in the 
top part of the right panel in Fig. [5] 

The situation changes in the case of a nucleus-nucleus collision: the bound nucleons in both nuclei participate in 
multiple interactions, therefore the scales of PDFs of all of them are modified. However, this modification goes beyond 
the simple shift Q 2 =>■ Q 2 + Q 2 A - Indeed, in an AB nuclear collision not only the two nucleons (one from A and one 
from B) participating in the hard reaction undergo multiple interactions, but also many other nucleons, the so called 
participants, experience multiple soft interactions. For this reason their parton distributions are boosted from the soft 
scale (i 2 up to the saturation scale fi 2 =>■ [i 2 + Q 2 mb)' wrnc h is usually much larger. Thus, the participant nucleons 
on both sides are boosted to a higher scale and get softer PDFs, with larger parton multiplicities at small x. This is 
illustrated on the cartoon in the bottom part of the right panel in Fig. El 



A. Reciprocity equations 

Thus, a participating nucleon simultaneously plays the roles of a beam and of a target. Its PDF is boosted to a 
higher scale due to multiple interactions it undergoes in another nucleus. As a target such a nucleon, being boosted 
to a higher scale Q\ => Q 2 , + Q 2 B , increases broadening of partons from another nucleus, since the factor C q (E, Q 2 ) 
Eq.@ rises. This leads to a mutual enhancement of the saturation scales in both nuclei. Indeed, multiple rescatterings 
of nucleons from the nucleus A on the boosted nucleons in B proceed with a larger cross section, so broadening, i.e. 
the saturation scale in B increases, Q 2 g Q 2 B > Q 2 s b- For this reason, the nucleon PDFs in A get boosted more. 
Then the partons from B experience even stronger multiple interactions with such double-boosted nucleons in A. 
This results in an additional boost of the saturation scale in A, then, as a result, in B, and so on. 

Such a multi-iteration mutual boosting of the saturation scales is illustrated pictorially in Fig. [6l where two rows 
of nucleons, Ta and Tb are displayed on horizontal and vertical axes. 
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FIG. 6: Collision of two one-dimensional rows of nucleons Ta and Tb displayed on horizontal and vertical axes. Multiple 
interactions of colliding gluons propagated through both nuclei, including additional multiple scatterings of gluons, which carry 
out the interactions, are shown as is described in text. 

To proceed further we present the function C(E,Q 2 ) in the form (45l. |4o|. 

C(E, Q 2 ) = y a s (Q 2 + Ql) xg N (x, Q 2 + Q 2 ), (11) 

where <7_/v(x, Q 2 ) is the gluon distribution function in the nucleon, and x = Im^Ejs is the fractional momentum of 
a gluon experiencing broadening. The scale shifting parameter Qo is introduced in order to have the correct infrared 
behavior: in the limit Q 2 — ?> Eq. (fTTJ) has to reproduce the known function C{E) depicted in the left panel of Fig. [5] 
This condition fixes the parameter Q 2 w 1.7 GeV 2 , which is nearly independent of energy. 
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Following Eq. (fTTj) the modified gluon saturation scales Q 2 in the collision of two rows of nucleons Ta and Tb can 
be found solving the reciprocity equations 44j |. 



37T 2 

1b Ob) = -rra s {Q 2 sA 



2 
3^ 
2 



<3o) xb9n{x b ,Q sA + Q ) T s; 
as(Q 2 sB + Qo) x A gN{xA, Q 2 sB + Ql) T A - 



(12) 



These equations take into account the modification of the properties of bound nucleons in each of the colliding nuclei 
due to multiple interactions in another nucleus and the following increase of the scale. 

Solutions of these equations [H| for xa = is = Qo/s for central collision of two identical nuclei at the energies of 
RHIC and LHC is plotted in the left panel of Fig. [7] as function of Ta, for the modified saturation scale (top) and 
for the ratio Ka = Qsa/QsA- For g(x, Q 2 ) we use the LO gluon distributions of the recent analysis MSTW2008 [43j . 
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FIG. 7: Left panel (top): the boosted values of the saturation momentum squared Q 2 3 a calculated at the energies of RHIC 
and LHC as function of nuclear thickness. Left panel (bottom): the boosting factor Ka = QIa/QsA as function of Ta at the 
energies of RHIC and LHC. Right panel: The boosting factors Ka(Ta) (solid curves) and Kb(Ta) (dashed curves) as function 
of nuclear thickness at ^/s = 5.5 TeV. Each pair of curves is marked by the rapidity for which it is calculated. 



We see that the saturation scale of heavy nuclei may be as large as about 10 GeV 2 at the LHC. In the right panel of 
Fig. [7] we plotted the ratios Ka(Ta) and Kb(Ta) with xa and xb calculated for charmonium production at different 
rapidities in the central collision of identical nuclei, A — B at ^fs — 5500 GeV [28j]. The saturation scales in A (solid 
curves) and B (dashed curves) are different, because the sign of rapidity is defined relative to the "beam" nucleus A. 



B. Boosted gluon shadowing 



As we have already mentioned, in terms of the Fock decomposition gluon shadowing corresponds to multiple 
interactions of higher Fock states, containing gluons. One may wonder of gluon shadowing in nuclear collisions 
factorize? 

Rf B (b,r) = Rf(r)Rf(b-f). (13) 

The answer is yes. The lifetime of gluonic fluctuation produced by a nucleon in the nucleus A may be sufficiently long 
only relative to the nucleus B, but is very short relative to the parent nucleus A. Therefore no gluonic fluctuations 
undergo double color filtering. In terms of Gribov inelastic shadowing this means that the diffractive excitation of 
the nucleons of A propagate through B independently of the excitations of B propagating through A. 
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One has to rely on a fully developed theoretical model for gluon shadowing to predict its 6-dependence. The effect 
of boosted saturation scale leads to a modification of the factor C(E) which is different for nuclei A and B, 

C(E) C A (E,T A )=K A (T A )C(E). (14) 

Solving equation ^ with such a boosted saturation scale one arrives at a modified nuclear ratio for gluons in nuclei 
A and B, R^ B \e, b), for which the factorized relation p^|) can be used (28j . 

Rf B (b,r) = Rf(r)R B (b-r). (15) 

A numerical example for the boosting effect on gluon shadowing is depicted by the solid curve in the right panel of 
Fig. [3] We see that the boosting effect considerably enhances gluon shadowing. 

V. FORWARD RAPIDITIES: CHALLENGING THE SUPPRESSION MECHANISMS 

A. Saturation effects 

The effects of saturation and shadowing are enhanced at small x, i.e. at higher energies. If the energy of collisions 
is fixed or restricted, one still can go to smaller x relying on the relation 

k 2 

xi x 2 = — , (16) 
s 

where 211,2 are not the values of Bjorkcn x for colliding partons, but are the fractional light-cone momenta of the 
produced high-fc^ parton relative to the beam and target respectively. Thus, one can reach smaller xi at fixed s, 
simply increasing x\ — > 1, i.e. moving to forward rapidities. Then the coherence phenomena (gluon shadowing, CGC) 
are expected to set up and suppress the particle production rate. However, as we multiply emphasized, these effects 
are rather weak, especially for the kinematics of RHIC experiments where the values of x 2 > 10 -3 , which can be 
accessed at forward rapidities, are not sufficiently small for gluon shadowing. 

Nevertheless, when a suppression was indeed observed in the BRAHMS experiment (34[, it was quite tempting to 
explain data by the effects of CGC. The model proposed in [47j contained few parameters fitted to the BRAHMS 
data. The results are depicted by the thick curve in Fig. [3] Notice that latest analysis [48[ also fits the absolute value 
of the ratio allowing different X-factors for nuclear and proton targets. Although agreement with data is reasonable, 
when data to be explained are fitted, the missed sources of suppression can be easily absorbed into the fit. One of 
them, the deficit of energy, is presented in the next section. 

B. Deficit of energy at forward rapidities 

Multiple interactions of the projectile hadron and its debris propagating through the nucleus should cause a dis- 
sipation of energy. This intuitive expectation is supported by consideration of the Fock state decomposition. The 
projectile hadron can be expanded over different states which are the fluctuations of this hadron. In the limit of 
infinite momentum frame those fluctuations live forever. One can probe the Fock state expansion by interaction with 
a target. The interaction modifies the weights of the Fock states, some interact stronger, some weaker. An example 
is the light-cone wave function of a transversely polarized photon (4lj |. In vacuum it is overwhelmed by qq Fock 
states with vanishingly small separation (this is why the normalization of the wave function is ultraviolet divergent). 
However, those small size fluctuations have a vanishingly small interaction cross section, and the photoabsorption 
cross section turns out to be finite. 

In each Fock component the hadron momentum is shared by the constituents, and the momentum distribution 
depends on their multiplicity: the more constituents are involved, the smaller is the mean energy per a constituent 
parton, i.e. the softer is the fractional energy distribution of a leading parton. So on a nuclear target the projectile 
parton distribution falls at x — > 1 steeper than on a proton. This is similar to the rescaling of the PDFs by the 
saturation scale we observed earlier (see Fig. [51 left panel), but a much stronger effect. 

In the case of a hard reaction on a nucleus, this softening of the projectile parton momentum distribution can be 
viewed as an effective loss of energy of the leading parton in the nuclear medium, because the initial state multiple 
interactions enhance the weight factors for higher Fock states in the projectile hadron. Those components with large 
number of constituents have a tough energy sharing, so the mean energy of the leading parton decreases compared to 
lower Fock states, which dominate the hard reaction on a proton target. Such a reduction of the mean energy of the 
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FIG. 8: Left panel: Nuclear modification factor RdAu of negative particles at r\ — 3.2. Data is from [34|. The solid line 
represents the main result of [43] for the h~ contribution. The meaning of other curves is explained in [43]. Right panel: 
Nuclear ratio, d — Au to pp, for negative particles as function of pr at pseudorapidity r\ = 3.2 [3J], and data for neutral pion 
production at r\ = 4 vs [3y] are depicted by round points and squares respectively. Solid and dashed curves correspond to 
calculations with the diquark size 0.3 fm and 0.2 fm respectively [43. |50|]. 



leading parton can be treated as an effective energy loss, which is proportional to the initial hadron energy. Indeed, 
the partons responsible for multiple Glauber collisions carry substantial fractions of the initial energy (formally, they 
should be included into the hadron-multipomeron vertex, rather than into a ladder [51|). Thus, the effective loss of 
energy is proportional to the initial energy. 

There is an important difference between this effect and energy loss of a single parton propagating through a 
medium and experiencing induced gluon radiation. In this case the mean fractional energy carried the radiated gluons 
vanishes with initial energy E as AE/E oc 1/E [521454] ]. This energy loss is independent of the parton energy. 

As far as multiple collisions suppress the cross section of leading particle production (heavy dilepton, high-p T 
hadron, heavy quarkonium, etc.), we assume that every collision brings in a suppression factor S(£) [43 ] . where 

£ = yjxl+xl; (17) 



2 PL 

XL = 7i ; 

xt = & (18) 

and pl.t are the longitudinal and transverse components of the momentum of the produced particle in the cm. frame. 
Notice that xl coincides with the Feynman variable, xl = xp. 

This factor S(£) should cause a strong (for heavy nuclei) suppression at £ — > 1, but some enhancement at small £ 
due to the feed down from higher £. This is because energy conservation does not lead to disappearance (absorption) 
of particles, but only their re-distribution in £. 

Since at £ — >• 1 the kinematics of an inelastic collision corresponds to no particle produced within the rapidity 
interval Ay ~ — ln(l — £), the suppression factor S(£) can be also treated as survival probability of a large rapidity 
gap, analogous to the Sudakov suppression for no gluon radiation. Assuming the Poisson distribution for the radiated 
gluons and using the gluon production rate estimated in [55| , the probability of such a rapidity gap was found [49j to 
be approximately, 

s(o « i - e (19) 
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This also goes along with the results of the dual parton model [56|, |57| . 

With such a suppression factor and applying the AGK cutting rules [51| with the Glauber weight factors, one 
achieves a parameter-free description of data depicted in the right panel of Fig. [5] With no adjustment the model 
also well describes the STAR data at 77 = 4 (36|. 

As an additional test, one can also look at the px spectra at different centralities (34J- The model agrees well with 
this data also, as is demonstrated in the left panel of Fig. |H1 
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FIG. 9: Left panel: Ratio of negative particle production in central (0-20%) and semi-central (30-50%) to peripheral (60-80%) 
d-Au collisions, shown by closed and open points respectively. The results of calculations [49J are plotted by solid and dashed 
curves. Right panel: Number of negative hadrons versus pt produced in pp collisions at y^s — 200 GeV and pseudorapidity 
n — 3.2. The results of the dipole model plotted by solid curve, are compared with BRAHMS data [34] ■ 



The last, but not least check is comparison with the px distribution of pions produced in pp collisions at this 
rapidity. The dipole model used in the calculations performed in [49( describes well both the absolute value and pr 
dependence of data [34| presented in the right panel of Fig. |H1 

Having two competing models, which are able to explain data at forward rapidities, one should look for specific 
reactions and kinematic domains where the models can be disentangled. Several tests are proposed below. 



1. Test #1: Down to smaller energies, other processes 



As we mentioned above, the main motivation for moving to forward rapidities is getting access to smaller Bjorken 
x. If smallness of x indeed is the reason for the observed suppression, this effect should disappear at lower energies, 
because x rises as oc 1/%/s- Having no other mechanism contributing to the suppression observed by BRAHMS, but 
only CGC, one should expect no suppression at forward rapidities at lower energies. 

On the other hand, the suppression caused by energy deficit scales in Feynman xf and should exist at any energy. 
Thus lowering the collision energy would be a sensitive test for the models. 

The NA49 experiment at SPS has performed measurements [HI similar to BRAHMS, but at much lower energy, 
where the value of xi is two orders of magnitude larger than in the BRAHMS data. The results depicted in the left 
panel of Fig. [TU] show that the effect of suppression at forward rapidities is still there. 

Moreover, it turns out that any reaction measured so far at forward rapidities exposes the same effect of an increased 
suppression. As an example, we present a collection of data for nuclear dependence (the exponent a fitted to data 
in the form of A a ) of the total yields of different species of hadrons in pA collisions as function of xf, in the right 
panel of Fig. 1101 The curve shows the result of the parameter-free calculation in (49j within the model described in 
Sect. EE 

Other reactions, like Drell-Yan process of heavy dilepton production, or charmonium production at forward ra- 
pidities also demonstrated a stronger suppression at large xf (see the data and references in [4!|). All these data 
have been taken at low energies of fixed target experiments, where no coherence effects can be expected. Apparently, 
another mechanism, which causes suppression at forward rapidities, is at work. It should be also important for the 
interpretation of RHIC data. 
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FIG. 10: Left panel: The nuclear ratio for pions produced in proton-lead collisions at Ei a i — 158 GeV as function of pr for two 
values of Feynman xf = 0.025 and 0.375 p^ - The curves present the results of the parameter free calculations [if| described 
in Sect. IV Bl Right panel: The exponent for the A a dependence fitted to data on total yield of different species of hadrons 
produced in pA collisions at different energies as function of xf ■ The curve present the results of the parameter free description 
of [U (see Sect.lVBj. 



2. Test #2: Towards large xt 



One can approach the kinematic limit £ — > 1 increasing either xp, or xt- In both cases the constraints imposed by 
energy conservation cause nuclear suppression. Therefore, we expect a suppression at large xt to be rather similar to 
what is observed for forward rapidities corresponding to large xl- This would be a rigorous test for the mechanisms 
of suppression, because no coherence is possible at large xt- 

The nuclear ratio R p a(pt) is known to expose nuclear enhancement at medium high px, the effect [sol ] named after 
James Cronin. The rather small magnitude ~ 10% of the Cronin enhancement for the dA to pp ratio, predicted 
in [37t | was confirmed by the later measurements of neutral pion production in the PHENIX experiment [39] , as is 
depicted by the dotted curves in Fig. [TTJ The expectation based on QCD factorization that the ratio must approach 
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FIG. 11: Nuclear attenuation factor RdAu{pT) as function of pr of ty° mesons produced in d-Au collisions at \fs = 200 GeV 
and 77 = 0. The solid and dashed curves represent the model predictions calculated with and without corrections Eq. (|19[) for 
energy conservation. The data [39( and calculations correspond to either minimum bias (left), or central (right) d-Au collisions. 

RdAuipr) — > 1 at large pr, is confirmed by the behavior of the dotted curves, with a small correction for the isotopic 
effects. 
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However, energy conservation becomes an issue with rising xx, and eventually causes a considerable suppression, 
as is shown by the solid curves in Fig. QTJ The same suppression factor Eq. (|19|) was used in this calculations, as at 
forward rapidities. The data do not contradict the predicted suppression at large px, even support it, especially in 
central collisions. Of course the observed suppression has nothing to do with the coherence effects. 

Notice, that the suppression at large px caused by energy deficit also contributes to the px dependence of hadron 
suppression in A A collisions. The observed flat px dependence of Raa{pt) would be rising, if the suppression caused 
by energy conservation were switched off. 

A similar test can be performed with direct photons produced with high px- This is even a cleaner probe, than 
hadron production, since photons have no final state interactions, and no convolution with the jet fragmentation 
function is required. Production of direct photons with large px approaching the kinematic limit also is subject to the 
energy sharing problem, so is universally suppressed by multiple interactions. Our predictions for high-p^ photons at 
the energies y/s — 200 and 62 GeV are depicted in Fig. [T2]in comparison with data from PHENIX experiment 60]. 
The energy range of RHIC is in between the two regimes of photon radiation, coherent and incoherent. To avoid the 
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FIG. 12: Ratio of cross sections of direct photon production in gold-gold to proton-proton collisions at y^s = 200 GeV (left 
panel) and at ^/s = 62 GeV (right panel) . Solid and dashed curves correspond to calculations done with and without the 
corrections for energy deficit, the factor Eq. (|19[1 suppressing multiple interactions, respectively. The data are from PHENIX 
experiment at RHIC [60| . 

technically complicated interpolation between the two regimes, we have done calculations for both regimes and found 
the difference to be quite small for px > 3 GeV. The solid and dashed curve are calculated for incoherent photon 
radiation, either including or excluding the corrections for energy limitations, respectively. Corrections for isotopic 
effects are included into these calculations. 

While nothing certain can be concluded from comparison with data at ^/s = 62 GeV, data at yfs = 200 GeV 
provide some evidence for the predicted suppression. In fact, the observed strong suppression of direct photons has 
been considered as a puzzle, since no energy loss or absorption in the produced hot medium is expected for photons. 
We explain suppression for particle production at large xl and large xx by the same mechanism. Of course the CGC 
scenario for forward rapidities [47j would have no effect on hadrons or photons produced with very large px- 

3. Test #3: The magnitude of gluon shadowing 

This test has been already done within the global fit [33| aimed at extraction of PDFs from data. The BRAHMS 
data [34| were included in the analysis, assuming that the observed suppression is related to gluon shadowing. The 
results depicted in the right panel of Fig. @] strongly violate the unitarity bound [27] ■ Apparently, the source of the 
trouble was the misinterpretation of the BRAHMS data. Moreover, even that strong shadowing was not sufficient to 
explain the STAR data at higher rapidity 77 = 4 [36| . 

VI. CRITICAL EXAMINATION OF THE ENERGY LOSS SCENARIO 



A parton propagating through a medium experiences multiple interactions, which induce gluon radiation taking 
away a part of the parton energy. Therefore, one can expect that the production leading hadrons, which carry the 
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main fraction of the jet energy will be suppressed compared to the same reaction in pp collisions [61j . Indeed a 
significant suppression, usually called jet quenching, was found at RHIC for all species of hadrons produced with large 
Pt in central collisions of heavy ions [62|, [63j ■ 

However, the energy loss scenario is based on the unjustified assumption that gluon radiation continues throughout 
the dense medium created in the collision, and the leading hadron (or a colorless state, called pre-hadron, which does 
not have any certain mass) is always produced outside the medium. 

This is quite a debatable issue, since there are solid theoretical and experimental arguments favoring the alternative 
scenario: the pre-hadron frequently is produced inside the medium and strongly attenuates. Indeed, it has been 
realized long time ago [64j that the production length for leading hadrons, may be rather short, 

1 >~WM {1 - Zh) - (20) 

Here E is the energy of the jet, dE/dl is the rate of vacuum energy energy loss, which is usually much larger than 
the medium induced one, especially for highly virtual partons. Here the suppression mechanism related to the deficit 
of energy is again at work. If the highly virtual parton intensively radiates, its energy steeply decreases (see next 
section) and can fall below the minimal energy required for production of the leading hadron of energy Eh = Zh E. 
Only creation of a colorless pre-hadron, which does not radiate gluons any more, can stop the dissipation of energy. 

In the string model the rate of energy loss —dE/dl = k ps 1 GeV/ fm is the string tension, known from the slope of 
the Regge trajectories, as well as from lattice calculations. The predicted z^-dependence of l p [64[ was confirmed by 
the study of time development of fragmentation within the Lund model [65| . 



A. Time dependence of vacuum radiation 

The color field of a quark originated from a hard reaction (high-px, DIS, e + e~, etc.) is stripped off, i.e. such a quark 
is lacking a color field up to transverse frequencies q <Q. Therefore the quark originated from such a hard process 
starts regenerating its field by radiating gluons, i.e., forming a jet. This can be described by means of an expansion 
of the initial "bare" quark over the Fock states containing a physical quark and different number of physical gluons 
with different momenta. Originally this is a coherent wave packet equivalent to a single bare quark \q). However, 
different components have different invariant masses and they start gaining relative phase shifts as function of time. 
As a result, the wave packet is losing coherence and gluons are radiated in accordance with their coherence times. 

This process lasts a long time proportional to the jet energy (Bsipj), since the radiation time (or length) depends 
on the gluon energy and transverse momentum k (relative to the jet axis), 



IE 2Ex(l .. , 

M^-m 2 k 2 +x 2 m 2 q ' V > 



Here x is the fractional light-cone momentum of the radiated gluon; m q is the quark mass; M 2 g = m 2 /(l — x) + 
k 2 /x(l — x) is the invariant mass squared of the quark and radiated gluon. 

One can trace how much energy is radiated over the path length L by the gluons which have lost coherence during 
this time interval [66i-[70j|, 

Q 2 l 

AE(L) =E [ dk 2 [ dxx-^-0(L - l c ), (22) 

A 2 

where Q ~ pr is the initial quark virtuality; the infra-red cutoff is fixed at A = 0.2 GeV. The radiation spectrum 
reads 



dn g _ 2a s (k 2 ) k 2 [l + (1 - x) 21 



dx dk 2 3tt x [k 2 + x 2 m 2 ] 2 



(23) 



where a s (k 2 ) is the running QCD coupling, which is regularized at low scale by replacement k 2 =>■ k 2 + k 2 , with 
k 2 = 0.5 GeV 2 . 

In the case of heavy quark the /c-distribution Eq. (|23[) peaks at k 2 x 2 m 2 , correspond ing to the polar angle (in 
the small angle approximation) 9 = k/xE = m q /E. This is known as the dead cone effect [Til. [72j|. 
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The step function in Eq. ([22} creates another dead cone: since the quark is lacking a gluon held, no gluon can be 
radiated unless its transverse momentum is sufficiently high, 



, 9 2Ex(l — x) o , 
k z > — - x 2 m q . 

Li 



(24) 



This bound is relaxing with the rise of L and reaches the magnitude k 2 ~ x 2 m 2 characterizing the heavy quark dead 
cone at 



E(l 



(25) 



We see that L q for beauty is an order of magnitude shorter than for charm, but linearly rises with the jet energy. 

The radiation of such a "naked" quark has own dead cone controlled by its virtuality Q 2 m 2 . This cone is much 
wider than the one related to the quark mass. There is no mass dependence of the radiation until the quark virtuality 
cools down to Q 2 =4> Q 2 {L) ~ m 2 . Therefore, the results of [72j for a reduced energy loss of heavy quarks should be 
applied with a precaution. At the early stage of hadronization, when Q 2 (L) ^> m 2 , all quarks radiate equally. 

The characteristic length L q may be rather long, since gluons are radiated mainly with small x. For instance, for 
a charm quark with E — px = lOGeV the sensitivity to the quark mass is restored at L > 1/x fm. Only at longer 
distances, L > L q , the dead cone related to the heavy quark mass sets up, and the heavy and light quarks start 
radiating differently. 

The numerical results demonstrating this behavior are depicted in the left panel of Fig. 1131 One can see that a 
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FIG. 13: Left panel: vacuum energy loss by light (m q — 0), charm (m c — 1.5 GeV) and bottom (mi, = 4.5 GeV) quarks with 
E — 15 GeV as function of path length. Right panel: Fractional energy loss in vacuum as function of distance covered by a 
light quark. The jet energies, E = pr, corresponding to the curves from top to bottom, are E = 5, 10, 15, and 20 GeV. 



substantial difference between radiation of energy by the charm and light quarks onsets at rather long distances, 
above 10 fm, while within several fermi the difference is insignificant. The 6-quark radiation is suppressed already at 
rather short distances. Moreover, it completely regenerates the color field already at a distance of the order of 1 fm 
and does not radiate any more. Of course, this &-quark still may have a medium induced radiation, which is very 
weak according to (72[. Notice that the interference between vacuum and induced radiations, which was found in [73| 
to be important, is absent because they occur on different time scales. 



B. Production length 



The fractional energy loss in vacuum by a light quark is depicted in the right panel of Fig. [T3] as function of path 
length. Right after the hard scattering the quark is highly virtual and the radiation is very intensive. As a result, 
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the quark virtuality is cooling down and the rate of energy loss decreases. So the quark dissipates most of energy at 
the early stage of hadronization, during the first 1 fm about 40% of the total radiated energy. A heavy quark does it 
even faster [70(. This means that a hadron with large fractional momentum Zh should be produced at a short time 
scale, otherwise too much energy will be radiated, and the rest will be insufficient for production of the hadron. This 
is what the equation (l20l) says. Notice that in a high-pT hadron production the large values of zt are favored by the 
steeply falling transverse momentum spectrum of quarks convoluted with the fragmentations function. 

In fact energy restrictions make the dissipation of energy for a given Zh somewhat slower than it follows from Eq. (|22[) 
and is depicted in F ig. | 13[ because radiation of gluons with fractional momenta larger than 1 — Zh is forbidden by 
energy conservation [6a, [69|. On the other hand, a ban for radiation of gluons with energy u > (1 — Zh)E in (|22l) . 
leads to the Sudakov type suppression factor, 

S(L, z h ) = exp [-<%(L, z h ))] , (26) 

where (n g (L, Zh)) is the mean number of nonradiated gluons during propagation over the distance L, 

l m aa= i 

(n g (L >Zh ))= j dl J da ^e(a + ^-l + , fc ). (27) 

l/Q (2B*)- 1 

Here a is the fractional light-cone momentum of a radiated gluon; l max = min{L, E/2X 2 }, and A is the soft cutoff 
for transverse momenta of gluons, fixed in [69( at A = 0.7 GeV. The step function in ([27| takes care of energy 
conservation. An example of a Sudakov factor calculated at E = Q = 20 GeV for different values of z is shown in the 
left panel of Fig. [HJ 
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FIG. 14: Left: Sudakov suppression caused by a ban for radiation of gluons with fractional energy higher than 1 — Zh- 
Calculations are done for a jet with E = 20 GeV. Right: The production length distribution dD(z)/dl p (in arbitrary units) for 
jet energies 6, 10, 16, 20 GeV and z h — 0.7. 



Combination of these two effects, vacuum energy loss and Sudakov suppression, leads to a rather short production 
length, which slightly varies with jet energy and virtuality. The resulting production length distribution for jets with 
maximal virtuality Q — pt/ %h = E is depicted in the right panel of Fig. [TJlfor different jet energies as was calculated 
in We see that the mean production length is rather short, few fermi, and is slowly decreasing with energy. We 
took into consideration so far only vacuum energy loss. Apparently, adding medium induced energy loss will only 
enhance the energy deficit and make the production length even shorter. 

The production length distribution calculated for light quarks is also valid for charm quarks, which have a similar 
vacuum radiation during first several fermi. However, a bottom quark, according to Fig. 1131 dissipates considerably 
less energy and its vacuum radiation ceases at the distance of about 1 fm, because the quark completely restores its 
color field. Of course according to confinement a colored quark even with a restored field cannot propagate freely. 
It keeps losing energy via nonperturbative mechanisms [69j, like in the mentioned above string (flux tube) model. 
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We can make a rough estimate of the production length for hadronization of a bottom quark, relying on Eq. (30] and 
assuming the pt distribution of the produced quarks do j dp\ oc p^ n , 

1 

( W =nf/*(l-z),^ = ^. (28) 
o 

This is quite a short distance, e.g. for E = 15 GeV and n = 6, it is only l p ss 2fm. 

Notice that one should not mix up the production time with the time scale evaluated in [74j . which is just the 
well known coherence time (e.g. compare Eq. (2) in 69] with Eq. (3) in 74 [ containing some misprints) , which is 
the time interval between the hadron production points where the production amplitudes interfere. This is not the 
time of duration of hadronization which we are interested in. If hadronization were lasting as long as the coherence 
time, energy conservation would be broken. Besides, a pre-hadron does not have any certain mass, since according to 
the uncertainty relation it takes time, called formation time, to resolve between the ground and excited states, which 
have certain masses. Therefore, one cannot evaluate the production time of a pre-hadron relying on the mass of the 
hadron. 

Thus, theoretically we see no justification for a long production length, while any reasonable evaluation leads to a 
rather short duration time for hadronization which ends up producing a leading hadron. Nevertheless, it would be 
more convincing to test the models comparing with data. 



C. Within or without? 



Strictly speaking this very question is not well defined. Indeed, in quantum mechanics one cannot always say 
with certainty whether the pre-hadron was produced within or without the medium. The corresponding amplitudes 
interfere like is illustrated in the left panel of Fig. [15J We employ here the Berger model for perturbative fragmentation 




FIG. 15: Left: Graphical representation for the interfering direct and conjugated amplitudes for production of a pre-hadron at 
points with coordinates Z2 and zz respectively. Right: Fractional cross sections <Tj/(<7i + 02 + 03) as function of zh calculated 
for lead at E = 10 GeV. 01,2,3 correspond to both coordinates 2:2,3 of pre-hadron production being either outside or inside the 
nucleus, or to their interference respectively [75| . 



[69j, l76l . |77| . The hard process occurs incoherently on a nucleon with longitudinal coordinate z\ . The radiated gluons 
decay into 5253 and q^q^ coherently in the two amplitudes. The colorless pre-hadrons (dipoles) 9291 and q^q'i created 
at zi and 2:3 respectively, are projected to the pion wave function in each of the two amplitudes. When both points 
Z2 and Z3 are located outside, or inside the nucleus, the corresponding parts of the cross section are labeled as a\ and 
<72 respectively. The interference term in the cross section, 173, corresponding to simultaneous production inside and 
outside, was found 100% important [75]. An example of relative contribution of all three term is shown in the right 
panel of Fig. 1X51 

Unfortunately, the phenomenology based on such a quantum-mechanical description has not been sufficiently de- 
veloped so far. Therefore, in what follows we treat the space-time development of in-medium hadronization in a 
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semi-classical way, and confront with data the two options: (i) depending on the process and kinematics the produc- 
tion time may be shorter or longer than the path length in the medium, correspondingly the hadronization ends inside 
or outside the medium; (ii) according to the basic assumption of the energy loss scenario the pre-hadron is always 
produced outside the medium. 



1. Test #1: Nuclear effects in SIDIS 



The shrinkage of the production length towards Zh = 1, Eq. (|2"U1) . and corresponding increase of nuclear suppression 
was predicted long time ago 1641 , but only the Hermes experiment performed the first measurements confirming such 
a behavior of nuclear ratios [78|. Moreover, the parameter-free quantitative predictions for the magnitude of nuclear 
attenuation, were well confirmed by later measurements (78|, as is demonstrated in the left panel of Fig. [TH] from [78[ 
for a nitrogen target. The main source of attenuation in the model (66j is absorption of the pre-hadron, which is a qq 
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FIG. 16: Left: HERMES data [78l ] for the multiplicity ratio as function of zu for all charged pions (open circles) and all charged 
hadrons including pions (closed squares). The full curve represents the prediction of the gluon-bremsstrahlung model [oq ]. 
Other curves correspond to different fits to the data. Right: The modified multiplicity ratios as function of Zh with different 
values of the jet transport parameter go {79t ] compared with the HERMES ^8fl| data for Ne, Kr and Xe targets. For clear 
presentation the modification factors for different targets are shifted vertically. 



dipole. Once produced inside the nucleus, it attenuates with a cross section, which is fluctuating during the formation 
of the hadron wave function. The fluctuations are summed up with the path integral technique. 

Because the production length Eq. (|20p linearly rises with energy (at a fixed virtuality), eventually the mean 
production length becomes so long, that any absorption corrections will be excluded, and the nuclear ratio should 
approach one. The energy dependence of nuclear attenuation predicted in (6(| was well confirmed by the HERMES 
measurements. At the energies of the EMC experiment v > 100 GeV no nuclear effects have been observed within the 
experimental errors. This is why the energy range of HERMES is most interesting and why this measurement were 
proposed in for HERMES. 

In the energy loss scenario, the whole effect of nuclear suppression at large Zh is attributed to medium induced 
energy loss, which modifies the splitting functions in the evolution equation and eventually the fragmentation function 
for hadronization in a medium. The authors admitted that the model failed to explain data at large Zh > 0.5 even by 
adjusting the transport coefficient qo, as is demonstrated in the right panel of Fig. [16] from [79| . 



2. Test #2: Heavy flavor production 

Heavy flavors have been predicted [72j to lose much less energy for in-medium induced radiation than light quarks, 
because of the dead-cone effect. Therefore they were expected to be less suppressed in heavy ion collisions. However 
in Sect. IVI Al we demonstrated that because of another, stronger dead cone, related to the very intensive vacuum 
radiation, charm quarks experience a reduced energy loss only on long path lengths, which are not important for 
quarks produced out of a dense medium. Otherwise, charm and light quarks radiate gluons with similar rates. This 
would be sufficient to explain the large suppression observed at RHIC high-pT electrons produced in central Au-Au 
collisions [8~il. |S2T | . However, recent measurements [83] found that an essential fraction of the electrons originates from 
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decays of 6-quarks. This fact leads to a serious problem for the energy loss scenario, which has failed so far to explain 
why beauty is suppressed as much as charm and light quarks [Hj|. 

On the other hand, we have just demonstrated in Eq. (|28p that the production length for beauty hadrons is very 
short, l p ~ 1 — 2fm, so the energy loss scenario is quite irrelevant. In this case absorption resulting from in-medium 
production is more important and is able to explain the observed strong suppression of 6-quarks as is discussed further 
in Sect. EED 



3. Test #3: Alternative probes for the dense medium 

Another problem of the energy loss scenario is too large density of the medium, which one needs to explain the 
observed jet quenching. For instance, the analysis performed in [85| led to the value of transport coefficient qo, which 
is more than order of magnitude larger than the conventional expectation [54| . This may indicate missed mechanisms 
of suppression, which were absorbed in the fit and may lead to wrong values of fitting parameters. A good test would 
be an alternative probe, which independently measures the same parameters. 

J production in nuclear collisions has been debated for a long time as a probe for the produced matter. It 
was predicted to dissolve at high temperatures as a result of Debye screening [8g|. However, no clear signal of this 
phenomenon has been observed so far, and eventually it was suggested [13] that it does not dissolve at all in the 
medium produced in heavy ion collisions in all previous measurements including RHIC. 

Nevertheless, a sizable suppression by break-up of J /'J' propagating through the created medium was found in 88]. 
The key point of this consideration is the relation between the dipole break-up cross section and pt broadening of a 
quark propagating through the medium [l5l [l6|. which provide direct access to the transport coefficient. 

In the left panel of Fig.[T7]we present the contributions of three different sources of nuclear effects for J production 
in Cu-Cu collisions, at y/s — 200 GeV. The curves from top to bottom show: (i) the contribution of J/^ absorption in 
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FIG. 17: Left: Nuclear ratio Raa for central (0-20%) copper-copper collisions as function of J /ty transverse momentum 
from the PHENIX (circles) and STAR [§(| (squares) measurements at yfs = 200 GeV. The dotted curve shows the FSI 
attenuation of J/\l/ in the produced dense medium. The dashed curve also includes the ISI effects, charm shadowing and 
absorption. The final solid curve is also corrected for the Cronin effect. Right: Nuclear ratio Raa for central copper-copper 
(full circles and squares, upper curve) and gold-gold (empty circles, bottom curve). 



the created medium, (ii) the nuclear suppression including initial state interactions (charm shadowing and break-up), 
and (iii) the final curve corrected for the Cronin effect. All the calculations are parameter free, except the upper 
(dotted) curve controlled by the transport coefficient, which was fitted and found q = 0.3 — 0.5 GeV 2 / fm. This value 
well agrees with the pQCD estimates |54|, but is significantly smaller than what results from the jet quenching analyses 
based on the energy loss scenario. The right panel of Fig. [T71 demonstrates the A dependence of J/^f suppression. 

It is worth mentioning that the predicted steep rise of Raa at pt > 5GeV is not reliable, since it is based on the 
parametrization of the pt dependence in pp collisions, which is extrapolated beyond the measured interval of pr- An 
alternative method of calculation of the Cronin effect, which was tried in [88[ led to the same results at px < 5 GeV, 
but a substantially weaker enhancement at larger px ■ 
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Notice, that the above analyses might miss some important mechanisms contributing to the observed nuclear 
suppression. One of them is the mutual boosting of the saturation scales in nuclei H . which results in an increased 
break-up cross section of J/ty. Another mechanism is related to the observation [9ll - l93j that the so called "cold 
nuclear matter" participating in the initial state interactions, is not really cold, as is discussed in (28[. Inclusion of 
these effects can only lead to a further decrease of the transport coefficient extracted from data. Therefore the above 
value of qo found from RHIC data, should be considered as an upper limit. 



VII. QUENCHING OF HIGH-p T HADRONS BY ABSORPTION IN A DENSE MEDIUM 

Saying that a pre-hadron produced inside the medium is absorbed, one should make it clear what it means, and 
why absorption suppresses the hadron production rate. Soft interactions cannot stop or absorb a high energy hadron 
or a parton, which only can change their color, while the longitudinal momentum remains unchanged. "Absorption" 
is a kind of jargon, which means that the interaction in the medium excludes or strongly suppresses the chance to 
detect the final state hadron in the given kinematic domain. This is not obvious for every interaction, e.g. in the case 
of break-up of a dipole. After a pre-hadron, i.e. a colorless dipole is produced it propagates without radiative energy 
loss, and the collision energy loss also is tiny, because the elastic cross section is very small. However, the dipole has 
a large color-exchange (inelastic) cross section. In the case of, say, charmonium dipole cc, a color exchange results 
in a break up and creation of two open charm hadrons. The chance to detect a charmonium becomes vanishingly 
small. However if the final hadron to be detected is a heavy-light quark dipole, say D-meson, originated from an 
initial charm quark, the break up is not so harmful, since the charm quark will always escape from the nucleus and 
produce a £>-meson with a high probability. Important, however, is to produce a leading Z?-meson with a certain 
and large fractional momentum Zh- This is actually why the production length shrinks at large Zh, Eq. (|20[) . as is 
depicted in Fig. Q31 The only way to stop the dissipation of energy is to produce a colorless pre-hadron. Once the 
latter breaks-up, the energy loss is again initiated and either a hadron with the given large z^ cannot be produced 
any more, or the Sudakov factor suppresses its production. 



A. How large is the pre-hadron? 

The pre-hadron propagating in a medium attenuates with a dipole cross section proportional to rj., where the mean 
transverse dipole separation tt is rising with time. A low energy dipole quickly expands to the hadronic size, but at 
high energies, Lorentz time dilation freezes the initial small size of the dipole for the time of propagation through the 
medium. So with rising energy of the dipole E, the medium becomes more transparent. 

The transverse expansion of a qq dipole moving can be described as, 

drr _ k T (t) 

dt a(l-a)E' [ ' 

where a is the fractional light-cone momentum of the quark. Applying the uncertainty relation kx ~ l/rr, we get, 

4(*)= n ^ , F +rl (30) 
a(l — a) E 

Where ro is the initial dipole separation. Such a behavior of the mean size squared can be also obtained within the 
more rigorous path integral technique [H|, [95[ for the early stage of expansion. 

The mean fractional momentum a can be estimated in the perturbative model for hadronization [7(1 W?\ , 

q^q + g* -> (qq) + q, (31) 

where the radiated virtual gluon decays to qq, and the q merges the parent quark creating a colorless qq dipole, which 
carries fraction of the initial quark momentum. This process is illustrated in Fig. 1151 Assuming equal momentum 
sharing in the g* — ¥ qq decay, one gets, 

a = 1. (32) 

Zh 

The expansion described by Eq. (|30p starts from the moment of the hard reaction which initiated the jet. The 
initial transverse extension of the color field in the high-p^ quark is very small, ro ~ 1/ptj and is quickly forgotten, 
since the second term in the right hand side of (|3U)) can be neglected. When the dipole pre-hadron is produced, its 
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mean initial size is given by Eq. (1301) at the distance from the hard interaction point I — l p . Then it keeps expanding 
in accordance with Eq. (|30l) . 

The size expansion at the early stage of jet development is so fast that the pre-hadron is produced with quite a 
large initial size, 

For example, for production length l p = 2fm, = 0.7 (the result is almost z^ independent for = 0.5 — 0.8), and 
jet energy E — pr/zh — 15 GeV, the dipole is produced with separation rx ~ 0.5 fm. This is a rather large size, which 
should cause a strong absorption in a dense medium. 



B. Attenuation of high-pr pions 

Let us consider a central, 6 = 0, collision of identical heavy nuclei with nuclear density PA( r ) — Pa®{Ra ~~ r )i 
where pa = 0.16 fm -3 . After the nuclei pass through each other the leave behind a quark-gluon medium, which is 
probed by produced high-px hadrons, as is illustrated in Fig. [18] Partons are produced in a hard collisions at impact 




FIG. 18: Left: Pictorial explanation of the space-time development of high-pr hadron production in a dense medium created 
in a central nuclear collision. Right: Pion suppression in central AA collisions (A ~ 200) at yfs = 200 GeV (solid) and 
yfs = 5500 GeV (dashed). Data are from the PHENIX experiment [g^]. 



parameter r with a rate proportional to Tf (t) = 4(R\ — t 2 )pa- Then the parton radiates gluons and propagates 
transversely through the expanding medium, whose density is diluting inversely with time. Besides vacuum radiation, 
the parton experiences multiple interactions in the medium which induce an additional radiation. At some point 
inside the medium, or outside, the parton picks up an anti-parton and produces a dipole, pre-hadron, which does not 
radiate any more. 

Integrating over the location of the hard collision and over t he p roduction time l p weighted with the distribution 
depicted in Fig. Q31 we arrive at the nuclear suppression factor [96f |97j , 



Raa = ^ 



R 2 



(34) 



Here the numerical factors a, /3 are independent of pt and the medium properties and are not important for further 
evaluations. The effective mean free path L has the form, 

L " = 8piR A X> (35) 

where the factor X besides other kinematic factors also includes the unknown initial density of the medium, which 
should be considered as a fitting parameter adjusted to data on R AA . However, there is a unique possibility to 
predict R AA without knowing the medium density. Assume that the medium density is so high, and the absorption of 
pre-hadron dipoles is so strong that their mean free path is very short, L <^ (l p ). In this case the survived pre- hadrons 
should be born either near the surface of the medium, or outside. Then, the further details of the medium features 
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are not important, since the jet can be produced only inside the outer slice of the medium of thickness (l p ), and this 
geometry defines the value of Raa- Indeed, in this case the last two terms in (I35[) can be neglected, and Raa much 
simplifies, 



This regime of hadron production assuming a very short mean free path L, looks quite plausible. As was demonstrated 
above, the starting transverse separation of the produced pre-hadron dipole given by Eq. (j3"3")l . is quite large and keeps 
expanding, what leads to a strong absorption in a dense medium. Indeed, the comparison of Eq. (|36p shown by the 
solid curve in the right panel of Fig.[T51 with RHIC data for pion production [62J demonstrates good agreement, both 
in magnitude and px dependence. This success encourages to believe that the high-density regime indeed occurs in 
the central collisions at RHIC. In this case it should be valid at the energies of LHC as well. Although a further 
increase of the medium density expected for LHC, should not affect Raa, Eq. (f31)|). the jet structure changes. While 
at RHIC the high pt pion production is dominated by valence quarks, at LHC gluons are expected to be the main 
source of jets. In this case the rate of perturbative vacuum radiation should be larger by the Casimir factor 9/4, and 
correspondingly (l p ) shorter. With the production length distribution calculated for gluon jets (9(| H3| the prediction 
for Raa at the energies of LHC is shown in Fig.[T5]by the dashed curve. Notice that as far as the medium is sufficiently 
dense for the regime of strong absorption, the dependence on collision energy should be very mild, most important is 
whether the hadron originates from quark or gluon jets. 

The validity of the assumption that the mean free path Eq. fl35[) is short depends not only on the medium density 
embedded in the parameter X, but also oupx- Because L is rising with px, eventually the approximation Eq. (|36[) will 
break down and Raa will start rising tending to restore the expectation based on QCD factorization, Raa 1- This 
hardly may happen at RHIC, because of the predicted suppression at large xt caused by the kinematic restrictions 
discussed in Sect. IV B 21 Nevertheless this should be expected for LHC. The position of the turn-over point in the pr 
scale should be sensitive to the medium density and deserves further study. 

Notice that in the case of production of two high-p^ hadrons absorption does not double, as one could naively 
expect. Raa will be similar to the case of single hadron detection. Indeed, in both cases the jet from which the one or 
two detected hadrons originate, must be within the outer medium slice of thickness equal to the shortest production 
length. Since l p does not vary much with the jet energy (see Fig. [T4")) , it does not make difference whether one, or two 
correlated hadrons are detected. 



As was already mentioned, a high-p^ heavy quark always escapes from the medium and produces an open flavor 
hadron. Such a process should have no suppression. Therefore a break-up of a light-heavy dipole in a medium does 
not lead to a suppression, unless the fractional momentum Zh of the detected hadron is fixed at a large value. In 
such a case break up of the dipole ignites continuation of vacuum energy loss, which slows down the quark to smaller 
values of z^. This is why a quark should stop radiating at a distance I ~ l p and produce a colorless dipole, which 
then survives through the medium. 

At first glance, since heavy quarks radiate less because of the dead-cone effect, their production length (l p ) should 
be longer compared with light quarks. However, as was pointed out in Sect. IVI Bl a charm and light quarks radiate 
similarly at short distances I ~ (l p ), therefore they should have similar production lengths. 

The main puzzle for the energy loss scenario, which remains to be a challenge, is the same strong suppression ob- 
served for bottom quarks. Although 6-quarks radiate much less than light ones [72|], according to the time-dependence 
of energy loss shown in the left panel of Fig. 1131 this occurs mainly because of prompt regeneration by 6-quarks of 
their color field. After that the 6-quark stops radiating, but forms a color flux tube [69| . which becomes the main 
source of energy loss. This leads to a quite short production length Eq. (|28p. 

Interesting that the produced heavy-light, c — q or b — q dipoles expand their sizes faster than a light qq dipole. 
This happens because of a very asymmetric sharing of the longitudinal momentum in such dipoles. Minimizing the 
energy denominator one gets the fractional momentum carried by the light quark, 



which indeed is very small, about 0.1 for charm and 0.03 for bottom. Then according to Eq. (l30l) a c — q dipole is 
expanding similar to qq, but a b — q dipole does it much faster. 

We conclude that hadronization of charm and bottom quarks ends up at a short distance l p with production of a 
colorless dipole. These dipoles are expanding similar, or even faster than a light qq dipole, therefore they are strongly 




(36) 



C. Heavy flavors 



m q 
a ~ — - 



(37) 
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absorbed by the dense medium. This justifies the application of the strong absorption scenario and Eq. (1361) to heavy 
charm and bottom quarks, and explains why both of them are strongly suppressed in A A collisions. 

VIII. SUMMARY 

This paper presents an attempt at a critical overview of the current status of our understanding of the dynamics 
of high-pT processes on nuclei. Here some of the important observations are highlighted. 

• The effects of coherence for gluon radiation are expected to be insignificant, since diffractive gluon radiation 
is suppressed in data. This means that gluon shadowing should be weak, what was confirmed by less biased 
global analyses of data for DIS on nuclei. Having a deeper insight at the contemporary global analyses, we 
concluded that the results for gluon nPDFs in some of them are not really constrained by data, but based on 
ad hoc assumptions, or otherwise rely on incorrect models for high-p^ hadronic reactions. 

• The same mechanisms which make gluon shadowing weak, also suppress the CGC effects. The saturation scale in 
nuclei, which is frequently overestimated, can be directly accessed by measuring px broadening of heavy flavors 
produced on nuclei. The dipole phenomenology fitted to DIS data correctly reproduces the data on broadening. 

• The bound nucleons in colliding nuclei considerably change their properties drifting to higher Fock components 
due to mutual multiple interactions. This boosts the saturation scale in AA compared with pA collisions as 
is described by the reciprocity equations. As a result, the nuclear medium becomes more opaque for colorless 
dipoles. A direct way to observe this effect in data is to compare the magnitudes of broadening observed in pA 
and AA collisions. 

• Trying to enhance the effects of coherence by decreasing Bjorken x, one should be cautious going to forward 
rapidities. Energy conservation becomes an issue towards the kinematic limits and may cause a strong suppres- 
sion. The same effect leads to a similar suppression of particles produced with large px (%t) hi pA and AA 
collisions, even if Xl is small. Indeed, data from RHIC for production of high-p^ pions and direct photons in 
dA and AA collisions provide evidence for such a suppression. 

• A highly virtual parton produced with high px i n a hard collision starts regenerating its color field and dissipating 
energy via gluon radiation. The energy loss is very intensive at the early stage, the parton radiates almost half 
of the radiated energy during the first 1 fm. Medium induced radiation speeds up this process, though usually 
it is a small correction. In order to respect energy conservation a leading hadron carrying a large fractional 
momentum Zh has to be produced at a short time scale after the hard reaction. This puts in doubt the main 
(unjustified) assumption of the energy loss scenario for jet quenching, that the hadronization is lasting a long 
time and ends up with production of a leading hadron always outside the medium. 

• The vacuum radiation of a highly virtual parton, which is lacking its color field, is subject to the dead-cone 
effect, making the radiation flavor independent during the initial stage of fragmentation. As a result, high-p^ 
charm and light quarks should be suppressed similarly in heavy ion collisions. Bottom quarks behave differently, 
they promptly restore the color field and stop radiating, on a distance of about 1 fm. The subsequent medium 
induced energy loss is too weak to explain the observed strong suppression of bottom produced in AA collisions. 
This remains a serious challenge for the energy loss scenario. 

• Neglecting the possibility of pre-hardon production and strong attenuation inside the medium, one should 
overestimate the medium density in order to reproduce the experimentally observed suppression of high-p^ 
hadrons. Indeed, this seems to be the case, the transport coefficient fitted to data on jet quenching within the 
energy loss scenario, substantially exceeds the conventional expectations. It is also much larger than the result 
of the alternative probe for the dense medium, suppression of J/*f? produced in AA collisions. 

• Once a pre-hadron is produced inside the medium, its subsequent fate depends on the initial dipole size and 
the speed of its expansion. Since the mean transverse momentum of the gluons radiated by a highly virtual 
parton steeply decreases with time, the initial transverse separation of the produced dipole rises as function of 
production length. For the mean production length we estimated that the produced pre-hadron dipole is rather 
large rx ~ 0.5 fm, and it keeps rising. Dipoles containing a heavy quark evolve even faster than light ones. As a 
result, the mean free path of such dipoles is very short, provided that the medium is substantially denser than 
the cold nuclear matter. 

This observation leads to a specific scenario for hadron production in a dense medium created in AA collisions. 
The hard interaction point must be located not deeper than at production length l p from the surface, otherwise 
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the produced pre-hadron will be promptly destroyed by the medium. As far as the model-dependent mean 
production length is known, this scenario predicts in a parameter- free way the suppression factor Raa, which 
well agrees with RHIC data. The suppression is similar for light and heavy quarks, including beauty. 



Acknowledgments 



B.K. is grateful to Lanny Ray for inviting to the Workshop and encouraging to write this paper. This work was 
supported in part by Fondecyt (Chile) grant 1090291, by DFG (Germany) grant PI182/3-1, by Conicyt-DFG grant 
No. 084-2009, by the Slovak Funding Agency, Grant 2/0092/10 and by Grants VZ MSMT 6840770039 and LC 07048 
(Ministry of Education of the Czech Republic) . 



0. V. Kancheli, Pisma Zh. Eksp. Teor. Fiz. 18, 465 (1973). 
L.D. Landau and I.Ya. Pomeranchuk, ZhETF 24, 505 (1953); 

L.D. Landau, I.Ya. Pomeranchuk, Doklady AN SSSR 92, 735 (1953); 
E.L. Feinberg and I.Ya. Pomeranchuk, Doklady AN SSSR 93, 439 (1953). 

L. D. McLerran and R. Venugopalan, Phys. Rev. D 49, 2233 (1994); Phys. Rev. D 49, 3352 (1994); Phys. Rev. D 50, 2225 
(1994). 

E. Iancu, A. Leonidov and L. D. McLerran, Nucl. Phys. A 692, 583 (2001) 

F. Gelis, T. Lappi and R. Venugopalan, Phys. Rev. D 78, 054019 (2008). 

L.V. Gribov, E .M. Levin and M.G. R yskin, Nucl. Phys. B188 (1981) 555; Phys. Rep. 100 (1983) 1. 

A. H. Mueller, |arXiv:hep-ph/9911289[ 

B. Z. Kopeliovich, L. I. Lapidus and A. B. Zamolodchikov, JETP Lett. 33, 595 (1981) [Pisma Zh. Eksp. Teor. Fiz. 33, 612 
(1981)]. 

B. Blaettel, G. Baym, L. L. Frankfurt and M. Strikman, Phys. Rev. Lett. 70, 896 (1993). 

B. Z. Kopeliovich, J. Raufeisen, A. V. Tarasov and M. B. Johnson, Phys. Rev. C 67, 014903 (2003) 

R. Baier, A. Kovner and U. A. Wiedemann, Phys. Rev. D 68, 054009 (2003). 

E772 Collaboration, E772 Collaboration, D.M. Aide et al., Phys. Rev. Lett. 66, 2285 (1991). 

M. B. Johnson, B. Z. Kopeliovich, M. J. Leitch, P. L. McGaughey, J. M. Moss, I. K. Potashnikova and I. Schmidt, Phys. 
Rev. C 75, 035206 (2007). 

P.L. McGaughey, J.M. Moss, and J.-Ch. Peng, NUCOLEX 99, Wako, Japan, 1999 < |hep-pli/9905447| ); Ann. Rev. Nucl. 
Part. Sci. 49, 217 (1999). 

M. B. Johnson, B. Z. Kopeliovich and A. V. Tarasov, Phys. Rev. C 63, 035203 (2001). 
J. Dolejsi, J. Hufner and B. Z. Kopeliovich, Phys. Lett. B 312, 235 (1993). 
B.Z. Kopeliovich, A. Schafer and A.V. Tarasov, Phys. Rev. D62 (2000) 054022. 

B. Z. Kopeliovich, J. Raufeisen, A. V. Tarasov and M. B. Johnson, Phys. Rev. C 67, 014903 (2003). 

Yu. M. Kazarinov, B. Z. Kopeliovich, L. I. Lapidus and I. K. Potashnikova, Sov. Phys. JETP 43, 598 (1976) [Zh. Eksp. 
Teor. Fiz. 70, 1152 (1976)]. 

A. Di Giacomo and H. Panagopoulos, Phys. Lett. B 285, 133 (1992). 

B. Z. Kopeliovich, I. K. Potashnikova, B. Povh and I. Schmidt, Phys. Rev. D 76, 094020 (2007). 
B. Z. Kopeliovich, J. Nemchik, A. Schafer and A. V. Tarasov, Phys. Rev. C 65, 035201 (2002). 

1. Balitsky, Nucl. Phys. B 463, 99 (1996). 

Y. V. Kovchegov, Phys. Rev. D 60, 034008 (1999). 

J. L. Albacete, N. Armesto, J. G. Milhano, C. A. Salgado and U. A. Wiedemann, Phys. Rev. D 71, 014003 (2005). 
B. Z. Kopeliovich, I. K. Potashnikova and I. Schmidt, Phys. Rev. C 81, 035204 (2010). 

in the right panel of Fig. [TS]B. Z. Kopeliovich, E. Levin, I. K. Potashnikova and I. Schmidt, Phys. Rev. C 79, 064906 
(2009). 

B. Z. Kopeliovich, I. K. Potashnikova, H. J. Pirner and I. Schmidt. larXiv: 1008.4272 1 [hep-ph] . 

A. Adare et al. [PHENIX Collaboration], Phys. Rev. C 77, 024912 (2008) [Erratum- ibid. C 79, 059901 (2009)]. 
D. de Florian and R. Sassot, Phys. Rev. D 69, 074028 (2004). 

K. J. Eskola, V. J. Kolhinen and P. V. Ruuskanen, Nucl. Phys. B 535, 351 (1998); 

K. J. Eskola, V. J. Kolhinen and C. A. Salgado, Eur. Phys. J. C 9, 61 (1999). 

K. J. Eskola, H. Paukkunen and C. A. Salgado, JHEP 0807, 102 (2008). 

I. Arsene et al. [BRAHMS Collaboration], Phys. Rev. Lett. 93, 242303 (2004). 

K. J. Eskola, H. Paukkunen and C. A. Salgado, JHEP 0904, 065 (2009). 

J. Adams et al. [STAR Collaboration], Phys. Rev. Lett. 97, 152302 (2006). 

B. Z. Kopeliovich, J. Nemchik, A. Schafer and A. V. Tarasov, Phys. Rev. Lett. 88, 232303 (2002). 
X. N. Wang, Phys. Rev. C 61 (2000) 064910. 

S. S. Adler et al. [PHENIX Collaboration], Phys. Rev. Lett. 98, 172302 (2007). 



27 



[40] M. Hirai, S. Kumano and T. H. Nagai, Phys. Rev. C 76, 065207 (2007). 
[41] B. Z. Kopeliovich, J. Raufeisen and A. V. Tarasov, Phys. Rev. C 62, 035204 (2000). 
[42] M. Arneodo et al. [New Muon Collaboration], Nucl. Phys. B 481 (1996) 3. 
[43] A. D. Martin, W. J. Stirling, R. S. Thorne and G. Watt, Eur. Phys J. C 64, 653 (2009). 
[44] B. Z. Kopeliovich, H. J. Pirner, I. K. Potashnikova and I. Schmidt. larXiv:1007.19l3l [hep-ph]. 
[45] L. Frankfurt, A. Radyushkin and M. Strikman, Phys. Rev. D 55 (1997) 98. 



[46] J. Bartels, K. J. Golec-Biernat and H. Kowalski, Phys. Rev. D 66, 014001 (2002) ;arXiv:he p-ph/0203258 . 
[47] D. Kharzeev, Y.V. Kovchegov and K. Tuchin, Phys. Lett. B599, 23 (2004). 
[48] J. L. Albacete and C. Marquet, Phys. Lett. B 687, 174 (2010). 

[49] B.Z. Kopeliovich, J. Nemchik, I.K. Potashnikova, M.B. Johnson,and I. Schmidt, Nucl. Phys. B146, 171 (2005). 

[50] J. Nemchik, V. Petracek, I. K. Potashnikova and M. Sumbera, Phys. Rev. C 78, 025213 (2008). 

[51] V. A. Abramovsky, V. N. Gribov and O. V. Kancheli, Yad. Fiz. 18 (1973) 595 [Sov. J. Nucl. Phys. 18 (1974) 308]. 

[52] F. Niedermayer, Phys. Rev. D 34 (1986) 3494. 

[53] S.J. Brodsky and P. Hoyer, Phys. Lett. 298B (1993) 165. 

[54] R. Baier et al., Nucl. Phys. B484 (1997) 265. 

[55] J.F. Gunion and G. Bertsch, Phys. Rev. D 25, 746 (1982). 

[56] A.B. Kaidalov, JETP Lett. 32, 474 (1980); Sov. J. Nucl. Phys. 33, 733 (1981); Phys. Lett. B 116, 459 (1982). 
[57] A. Capella et al., Phys. Rep. 236, 225 (1994). 

[58] B. Boimska, Ph.D. Dissertation (Warsaw 2004), CERN-THESIS-2004-035. 
[59] J.W. Cronin at al., Phys. Rev. D 11, 3105 (1975). 

[60] T. Isobe et al, [PHENIX Collaboration], J. Phys. G 34, S1015 (2007); T. Sakaguchi et al, [PHENIX Collaboration], Nucl. 

Phys. A 805, 355 (2008). 
[61] M. Gyulassy and M. Pliimer, Nucl. Phys. A 527, 641 (1991). 
[62] S. S. Adler et al. [PHENIX Collaboration], Phys. Rev. Lett. 91, 072301 (2003). 
[63] J. Adams et al. [STAR Collaboration], Phys. Rev. Lett. 91, 072304 (2003). 

[64] B.Z. Kopeliovich and F. Niedermayer, Sov. J. Nucl. Phys. 42 (1985) 504; Yad. Fiz. 42 (1985) 797. 
[65] A. Bialas and M. Gyulassy, Nucl. Phys. B 291, 793 (1987). 

[66] B.Z. Kopeliovich, J. Nemchik and E. Predazzi, Proceedings of the workshop on Future Physics at HERA, ed. by G. 

Ingelman, A. De Roeck and R. Klanner, DESY 1995/1996, v. 2, 1038; arXiv: |nucl-th/960"7036| 
[67] B. Z. Kopeliovich, J. Nemchik, E. Predazzi and A. Hayashigaki, Nucl. Phys. A 740, 211 (2004). 

B. Z. Kopeliovich, J. Nemchik and I. Schmidt, Nucl. Phys. A 782, 224 (2007). 
[69] B. Z. Kopeliovich, H. J. Pirner, I. K. Potashnikova and I. Schmidt, Phys. Lett. B 662, 117 (2008). 
[70] B. Z. Kopeliovich, I. K. Potashnikova, and I. Schmidt, arXiv: 1004.3920; to appear in Phys. Rev. C. 
[71] Y. L. Dokshitzer, V. A. Khoze and S. I. Troian, J. Phys. G 17, 1481 (1991). 
[72] Y. L. Dokshitzer and D. E. Kharzeev, Phys. Lett. B 519, 199 (2001). 
[73] N. Armesto, C. A. Salgado and U. A. Wiedemann, Phys. Rev. D 69, 114003 (2004). 
[74] C. Markert, R. Bellwied and I. Vitev, Phys. Lett. B 669, 92 (2008) [arXiv:0807.T509l [nucl-th]] . 

[75] B. Z. Kopeliovich, H. J. Pirner, I. K. Potashnikova, I. Schmidt, A. V. Tarasov and O. O. Voskresenskaya, Phys. Rev. C 

78, 055204 (2008). 
[76] E. Berger, Z. Phys. C4 (1980) 289. 

[77] B. Z. Kopeliovich, H. J. Pirner, I. K. Potashnikova, I. Schmidt and A. V. Tarasov, Phys. Rev. D 77, 054004 (2008). 



[78] A. Airapetian et al. [HERMES Collaboration], Eur. Phys. J. C 20, 479 (2001) arXiv:hep-ex/0012049 

[79] W. t. Deng and X. N. Wang, Phys. Rev. C 81, 024902 (2010) l arXiv:0910.3403l [hep-ph]]~ 

[80] A. Airapetian et al. [HERMES Collaboration], Nucl. Phys. B 780, 1 (2007) [arXiv: 0704 .32701 [hep-ex]] 

[81] A. Adare et al. [PHENIX Collaboration], Phys. Rev. Lett. 98, 172301 (2007). 

[82] B. I. Abelev et al. [STAR Collaboration], Phys. Rev. Lett. 98, 192301 (2007) 

[83] A. Adare et al. [PHENIX Collaboration], Phys. Rev. Lett. 103, 082002 (2009). 

[84] S. Wicks, W. Horowitz, M. Djordjevic and M. Gyulassy, Nucl. Phys. A 783, 493 (2007). 

[85] A. Adare et al. [PHENIX Collaboration], Phys. Rev. C 77, 064907 (2008). 

[86] T. Matsui and H. Satz, Phys. Lett. B178 (1986) 416. 

[87] F. Karsch, D. Kharzeev and H. Satz, Phys. Lett. B 637, 75 (2006). 

B. Z. Kopeliovich, I. K. Potashnikova and I. Schmidt, Phys. Rev. C 82, 024901 (2010). 
[89] A. Adare et al. [PHENIX Collaboration], Phys. Rev. Lett. 101, 122301 (2008). 
[90] B. I. Abelev et al. [STAR Collaboration], Phys. Rev. C 80, 041902 (2009). 
[91] J. Hufner and B. Z. Kopeliovich, Phys. Lett. B 445, 223 (1998). 
[92] J. Hufner, Y. B. He and B. Z. Kopeliovich, Eur. Phys. J. A 7, 239 (2000) 
[93] J. Hufner, B. Z. Kopeliovich and A. Polleri, Eur. Phys. J. A 11, 457 (2001) 
[94] B. Z. Kopeliovich and B. G. Zakharov, Phys. Rev. D 44, 3466 (1991). 
[95] B. Z. Kopeliovich, A. V. Tarasov and A. Schafer, Phys. Rev. C 59, 1609 (1999). 

[96] N. Armesto et al, J. Phys. G 35, 054001 (2008). 

[97] B. Z. Kopeliovich, I. K. Potashnikova and I. Schmidt. larXiv:0707.43"02l [nucl-th]. 



